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Abstract 

The widespread use of metallic structures in space technology brings risk of degradation which occurs under space 
conditions. New types of materials dedicated for space applications, that have been developed in the last decade, are in 
majority not well tested for different space mission scenarios. Very little is known how material degradation may affect 
the stability and functionality of space vehicles and devices during long term space missions. 

Our aim is to predict how the solar wind and electromagnetic radiation degrade metallic structures. Therefore both 
experimental and theoretical studies of material degradation under space conditions have been performed. The studies 
are accomplished at German Aerospace Center (DLR) in Bremen (Germany) and University of Zielona Gora (Poland). 

The paper presents the results of the theoretical part of those studies. It is proposed that metal bubbles filled with 
Hydrogen molecular gas, resulting from recombination of the metal free electrons and the solar protons, are formed on 
the irradiated surfaces. A thermodynamic model of bubble formation has been developed. We study the creation process 
of H 2 -bubbles as function of, inter alia, the metal temperature, proton dose and energy. Our model has been verified by 
irradiation experiments completed at the DLR facility in Bremen. 

Consequences of the bubble formation are changes of the physical and thermo-optical properties of such degraded 
metals. We show that a high surface density of bubbles (up to 10 8 cm -2 ) with a typical bubble diameter of ~ 0.4pm will 
cause a significant increase of the metallic surface roughness. This may have serious consequences to any space mission. 

Changes in the thermo-optical properties of metallic foils are especially important for the solar sail propulsion 
technology because its efficiency depends on the effective momentum transfer from the solar photons onto the sail 
structure. This transfer is proportional to the reflectivity of a sail. Therefore, the propulsion abilities of sail material 
will be affected by the growing population of the molecular Hydrogen bubbles on metallic foil surfaces. 
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1. Introduction 

Metallic structures are commonly used in space tech¬ 
nology. They build skeletons of spacecrafts, they protect 
satellites’ interiors from rapid temperature changes (MLI 
blankets), or they are used as highly reflecting mirrors of 
optical space telescopes. Nowadays, metals are also used 
as thin layers on polyimide-type foils which have a broad 
usage in the solar sail technology. 

A failure of a space mission may be the result of a 
change of metallic structure properties caused by the en- 
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vironmental effects. Therefore, all of the materials planned 
for space applications have to be evaluated for their behav¬ 
ior under particle and electromagnetic radiation (ECSS-Q- 
ST-70-06C, 2008; ASTM, 2010). It is known from many 
of these evaluation tests that particle and electromagnetic 
radiation can significantly degrade materials (see e.g. Lura 
et al., 1997; Heltzel et al., 2009; Sharma & Sridhara, 2012). 

Many kinds of metals have found their application in 
space industry, e.g. Aluminum, Copper, Nickel, Tita¬ 
nium, steels, and others (ECSS-Q-70-71A, 2004; ECSS-Q- 
70-71C, 2014). For instance Aluminum is one of the basic 
building materials of existing spacecrafts and is compo¬ 
nent of many subsystems. Copper is used in electrical, 
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electronic and in general engineering applications. Nickel 
has its application e.g. in heating elements. Titanium is 
chosen in space applications for its mechanical-, temper¬ 
ature - properties, and chemical resistance (ECSS-Q-70- 
71A, 2004; ECSS-Q-70-71C, 2014). 

However, it is very clearly stated that the radiation at the 
level existing in space does not modify the properties of 
metals (ECSS-Q-70-71A, 2004). 

The here presented paper proves that the thin metal¬ 
lic foils are especially sensitive to the ion irradiation. The 
free electrons within the metals can in a well defined energy 
range recombine with the solar wind protons into neutral 
Hydrogen atoms. Recombination processes and their con¬ 
sequences have direct effects onto the foils’ physical and 
thermo-optical properties (Sznajder & Geppert, 2014). 

The result of a recombination of metal free electrons 
and solar protons is a formation of bubbles filled with Hy¬ 
drogen molecular gas (see e.g. Milacek et al., 1968; Szna¬ 
jder & Geppert, 2014). Bubble formation is one of the four 
degradation mechanisms caused by Hydrogen (referred to 
as embrittlement): formation of a hydride phase, enhanced 
local plasticity, grain boundary weakening and bubble for¬ 
mation (Myers et al., 1992; Lu & Kaxiras, 2005). 

When H 2 -bubbles are formed on irradiated metallic 
surfaces, their reflectivity decreases with increasing sur¬ 
face density of the bubbles (Sznajder & Geppert, 2014). 

A reflectivity is a key parameter for many thin metallic 
foil applications e.g. in the solar-sail propulsion technol¬ 
ogy- 

Since the reflectivity is directly proportional to the mo¬ 
mentum transfer from solar photons to the sail’s material, 
its decrease will reduce the propulsion performance of a 
sail-craft just linearly proportional to the reduction of the 
reflectivity (Sznajder & Geppert, 2014). 

The future solar sail missions will be realized in the in¬ 
terplanetary medium (see e.g. Geppert et al., 2011; Kawaguchi 
J., 2014; Macdonald M. et al., 2014). Unfortunately, the 
real degradation behavior of metallic samples is to a great 
extent unknown. Hence, detailed studies, both theoreti¬ 
cal and experimental, are performed at German Aerospace 
Center (DLR) in Bremen, Germany with a cooperation of 
University of Ziclona Gora, Poland. 

The here presented thermodynamic model simulates 
the growth of H 2 -bubbles. The model input parameters 
are: the energy and flux of solar protons, type, and the 
temperature of the irradiated metal. The diffusivity of 
H in the metal lattice was taken into account, as well as 
back scattering effect (BS) of the solar protons irradiating 
the target. The model output is the velocity of bubble 
radius growth, the maximum possible bubble radius, and, 
for a given bubble density and average bubble radius, the 
reduction factor of the reflectivity with respect to its ideal 
value. 

The paper is organized as follows. In Section 2, the 
typical proton flux spectra at the distance of 1.0 AU from 
the Sun are presented. The value of the proton flux can 
be transformed to any distance orbit from the Sun, see 


Eq. 2. Four recombination processes of the metal free 
electrons and the incident protons are presented: Auger 
resonant-, Oppenheimer-Brinkman-Kramers (OBK), and 
Radiative Electron Capture (REC) - process. In Section 3 
general principles and conditions for the H 2 -bubble forma¬ 
tion are given. Next, the thermodynamic model of bubble 
formation is introduced. Also a simple model of specular 
reflectivity reduction due to growing population of bub¬ 
bles is presented. In Section 4 experimental results as well 
as validation of the thermodynamic model are presented. 
The reflectivity of metallic surfaces which are populated 
by a given surface density of bubbles is studied as well. 
Finally, in Section 5, the conclusions are drawn. 

2. Recombination of protons into neutral Hydro¬ 
gen atoms 

Devices, while operating in the interplanetary space, 
are exposed to solar wind and electromagnetic radiation. 
The solar wind, as the Sun’s corona, is essentially made up 
of electrons and protons plus a small proportion of heavier 
ions, and it carries a magnetic field. Particles and fields 
are intimately coupled in plasmas (Meyer, 2007, Ch.l). 

Extraterrestrial Sun observatories measure a few key 
solar wind parameters, e.g.: components of proton and 
electron velocity, their mean number density as well as 
components of the magnetic field. 

Fig. 1 shows solar proton at 1 AU distance from the 
Sun for its average activity. Proton fluxes are calculated 
by use of the data collected by the SOHO (since 1995) and 
the ACE (since 1997) satellites. The OMERE database is 
also considered. 


Bid 


ACE database (p+) 

SOHO database (p+) 

OMERE database - GEO Orbital Average FI 


E [keV] 


Figure 1: Flux of solar protons as a function of energy. Data are 
taken from the SOHO, ACE, and OMERE database. 


When a probe is irradiated in space, it collects inci¬ 
dent ions from a wide energy range. The range depends on 
the type and the thickness of an irradiated material. The 
thiner the target material, the less ions stuck in it. There¬ 
fore, there must exist a critical energy of incident ions ( E c ) 
above which they pass through the material. Hence, the 
integrated proton flux over the energies is: 
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where E m j n is the ion’s lowest energy recorded by the satel¬ 
lite’s detector system. The Ie values are presented in the 
Table 1. To calculate the fluxes the ACE database was 
used. 

Obviously, the E c depends on materials’ type and thick¬ 
ness. The 7.5 /iin Upilex — S® foil covered on both sides 
with 100 nm A1 layer has been used to perform the first 
proton irradiation tests at DLR. For that coating thickness 
all of the protons with energy lower than 8 keV will stuck 
in the Aluminum layer. Therefore, the integrated proton 
flux is 1.15 x 10 13 [p + cm -2 s -1 ]. 

Table 1: Integrated proton fluxes over the energies for 1 AU distance 
orbit from the Sun. 


E c [keV] 

Ie x 10 13 [p + cm 1 2 s x ] 

1.0 

0.44 

1.5 

0.68 

2.0 

0.91 

2.5 

1.06 

3.0 

1.12 

4.0 

1.14 

5.0 

1.15 

9.0 

1.15 


To estimate the flux of solar protons Id at distance d 
from the Sun, the following relation can be used: 

4tt( 1AU) 2 x I E = 4?rd 2 x I d . (2) 

Incident protons, while penetrating the metallic tar¬ 
get, recombine with its free electrons to neutral Hydrogen 
atoms. There are four recombination processes of ions into 
neutral atoms: 

1. Auger process 

In the Auger processes, an electron is captured by the 
incident ion to a bound state (here a proton) forming 
neutral Hydrogen. In order to conserve energy an 
Auger-electron or photon is released (see e.g. Guinea 
et al., 1982; Sols & Flores, 1984, 1988; Echenique et 
al. , 1990; Penalba et al., 1990; Roesler & Garcia de 
Abejo, 1996; Pauly et al., 2002). 

2. Resonant process 

The resonant recombination proceeds when the inci¬ 
dent ion is neutralized by an electron which is tun¬ 
neled to a metastable state (Hagstrum, 1954). The 
inverse process is also possible. An electron which 
is in a metastable state with respect to the metallic 
ion can populate one of the free electron states of 
the metal only if it becomes free (the Pauli exclusion 
principle). 

The effect comes from the crystal structure itself. 
The resonant processes are due to the potential seen 


by the moving ion i.e. they are described in a frame 
where the incident ion is at rest (Sols & Flores, 1988). 
From the point of view of the ion, there appear a 
moving periodic potential which gives rise to transi¬ 
tions between bound states of the ion-electron pair 
(composite) and free electron states (Echenique et 
al. , 1990). 

3. Oppenheimer-Brinkman-Kramers (OBK) process, 
The OBK process is a capture process, where an in¬ 
ner or outer shell electron of a target atom is trans¬ 
ferred to the moving ion (Sols & Flores, 1988). In 
the literature there are many physical approaches 
(see e.g. Chew, 1950; Bransden & Cheshire, 1963; 
Lapicki & Losonsky, 1977; Lin et al., 1978; Belkic et 
al., 1979; Ford et al., 1981; Alston, 1983; Miraglia, 
1984; Ghosh et al., 1987; Gravielle & Miraglia, 1988; 
Decker & Eichler, 1989; Datta et al., 1990). Different 
results may be obtained depending on the approxi¬ 
mation applied to the wave functions and the energy 
levels involved in the process (Echenique et al. , 

1990) . 

In the OBK process the outer-shell electrons of metal 
ions experience a strong Coulomb field of the inci¬ 
dent ion. The wave function of the electrons is dis¬ 
torted (Kuang I, 1991; Kuang II, 1991). For the 
inner-shell capture, the screening effect of the outer- 
shell electrons of the metallic ions reduces the cap¬ 
ture probability of an electron by the incident par¬ 
ticle (Winter & Lin, 1974; Banyard, 1977; Kuang I, 

1991) . 

4. Radiative Electron Capture (REC) process. 

In the REC process an electron is transferred from 
the target atom to the incident ion with the simulta¬ 
neous emission of a photon, see (Eichler & Stoehlker, 
2007) or (Eichler, 2005, Ch. 10, p. 151). It occurs 
at ion energies exceeding 150 MeV (for Aluminum 
as target material) (Raisbeck & Yiou, 1971). 

Since the solar wind consists mainly of low (< 100 keV) 
energetic protons, only the first three processes are respon¬ 
sible for recombination. 

The efficiency of the recombination processes is deter¬ 
mined by their cross sections E. The cross section of each 
recombination process is calculated by use of the concept 
of so-called transition amplitude A t . It determines the 
probability, P = A t , for a transition from an eigenstate 
Q to Q' (Penrose & Isham, 1986, Ch. 15, p. 227). Here 
P stands for probability of a recombination event i.e. that 
an electron is bound to an incident ion. Aj? is the complex 
conjugate of A, Q is an eigenstate of the incident ion and 
an electron before recombination takes place, while Q' is 
an eigenstate of the ion-electron composite. Hence, the 
cross section of recombination event can be written as: 

pOO 

E ~ / A* A t b db, (3) 

Jo 

where b is the impact parameter. 


3 





Fig. 2 shows the different cross sections for capture 
processes when Aluminum is irradiated with protons (H + /A1) 
as a function of proton kinetic energy given in keV (Sols 
& Flores, 1984). One can see that the Auger process £ a is 
the dominant one. The resonant process Sr has negligible 
contribution to the total cross section. The OBK process 
£obk (in the literature it is often called the shell process) 
gives the main contribution to the total cross section at 
some MeV, depending on the target material (Raisbeck & 
Yiou, 1971). 
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Figure 2: Cross sections for capture processes of H+ ion: Auger (Ea), 
resonant (Sr) and OBK (Sqbk) process. The dashed line represents 
the total cross section S tota i of all processes (Sols & Flores, 1984). 


3. Formation of molecular Hydrogen bubbles 

H 2 -bubbles occur as an irradiation damage. They change 
the physical properties of the irradiated surface and in¬ 
crease the erosion rate (Astrelin et al., 2010). H 2 -bubbles 
are metal pockets filled with Hydrogen molecular gas. The 
tendency to form bubbles depends on the proton energy, 
integrated proton flux (dose), temperature of the target, 
crystallographic orientation of the irradiated surface as 
well as on impurities and defects in the sample (Daniels, 
1970). It is known from terrestrial laboratory experiments 
that for Aluminum the minimum dose of protons above 
which the process occurs is ~ 10 16 H + cm~ 2 (e.g. Milacek 
et al., 1968). The temperature range in which bubbles 
were observed is between 288 and 573 K (e.g. Milacek et 
al., 1968; Daniels, 1970). 

The procedure that was used to estimate the critical 
temperature (573 K) above which the process of bubble 
formation was stopped due to the bubble cracking mecha¬ 
nism was as follows. The Aluminum target was irradiated 
by a flux of protons at room temperature. When irradi¬ 
ation of the sample was stopped, the probe was heated 
up to higher temperatures. A significant increase of both, 
the surface density and sizes of the bubbles has been ob¬ 
served until the critical temperature was reached. That 


procedure, used by the authors, allows to capture more 
Hydrogen by the vacancies since during the irradiation, 
and at room temperature the vacancies will collect more 
Hydrogen than at elevated temperatures. The vacancy is 
a missing ion or point defect in metal lattice (Damask & 
Dienes, 1971). Also a diffusion of Hydrogen in Aluminum 
at room temperature is much lower than at temperatures 
reaching ~ 570 K (Linderoth, 1988). In space a probe is 
bombarded by the solar protons at the temperature which 
is related to its orbit. Therefore, the procedure presented 
by (Milacek et al., 1968; Daniels, 1970) does not match the 
bubble formation mechanism under real space conditions. 

Hydrogen atoms are much smaller than metal ions, but 
they can introduce strain in a metal lattice when being ab¬ 
sorbed as interstitial ions (Matzger et al., 1976; Thomas 
& Drotning, 1983; Ren et al., 2008). They can also change 
the electronic structure of near neighbor metal ions (Ren 
et al., 2008). That causes an increase of the lattice energy. 
It may be decreased by the aggregation of the interstitial 
Hydrogen atoms into Hydrogen atom clusters, and then 
molecular Hydrogen bubbles (Ren et al., 2008). Hydro¬ 
gen could not agglomerate into ^-clusters without the 
presence of vacancies. For instance a single vacancy in 
Aluminum can trap up to twelve H atoms. For compar¬ 
ison, a vacancy in Iron can trap only up to six H atoms 
(Lu & Kaxiras, 2005). Hydrogen atoms which do not con¬ 
stitute the H 2 molecules and which are not trapped by the 
vacancies, diffuse through the metal lattice. 

Molecular Hydrogen bubbles were observed also on dif¬ 
ferent materials than Aluminum. Copper, Tungsten, Pal¬ 
ladium and Iron were investigated (e.g. Astrelin et al., 
2010). Bubbles are not forming on Tantalum and Vana¬ 
dium. These metals are well known as blistering-resistant 
materials (Astrelin et al., 2010). However, they are not 
suitable for space applications where the surface reflec¬ 
tivity plays a crucial role, e.g. in the solar sail propulsion 
technology since their reflectance is ~ 50% lower than that 
of Aluminum (Polyanskiy, 2014). 

3.1. Formation of bubbles under space conditions 

Growth of molecular Hydrogen bubbles will be possible 
in the interplanetary space if the criterion of the minimum 
dose of protons is fulfilled. The temperature of the sample 
has to be high enough to start the bubble formation, but 
not too high to lose Hydrogen much too rapidly due to the 
high diffusivity of Hydrogen in metals. 

Under the simplifying assumption that the Sun gener¬ 
ates only mono-energetic 5 keV protons, the criterion of 
minimum dose of protons will be fulfilled after 116 days 
for 1 AU distance orbit from the Sun. Obviously, taking 
into account proton fluxes from the whole energy range, 
the criterion will be fulfilled much earlier. 

The temperature of a foil placed in a given distance d 
from the Sun can be calculated from the balance of heating 
and cooling by: 
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Figure 3: Temperature of the Upilex — S® foil covered on both sides 
with 100 nm vacuum deposited Aluminum layer as a function of the 
distance from the Sun. The red area represents temperature range 
in which the bubble formation was reported in the literature. The 
dark-red area is the temperature range in which the formation has 
been confirmed by studies performed at DLR. 
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Here, A a is the area of the sample which absorbs the elec¬ 
tromagnetic radiation, while A e is the area which emits 
the heat by radiation. Hence, the ratio ^ equals 0.5. 
(Jsb is the Stefan-Boltzmann constant, SC states for So¬ 
lar Constant. The thermo-optical parameters have been 
provided by the manufacturer of the Upilex — S® foil, the 
UBE company. Solar absorptance as and normal emit- 
tance et are 0.093 and 0.017, respectively. The foil tem¬ 
perature as a function of distance from the Sun is repre¬ 
sented by solid line in Fig. 3. Note, that the heat released 
by stopped protons is negligible small in comparison to the 
Sun’s input. The red area (570 - 300 K) is the temperature 
range in which the bubble formation has been confirmed 
by the terrestrial laboratory experiments. Unfortunately, 
commonly used experimental procedures to estimate the 
maximum temperature at which the bubble formation is 
stopped, are not suitable for the real space conditions. The 
dark-red area represents temperatures at which the bub¬ 
ble formation has been confirmed by the first experimental 
findings performed at DLR. Some of the results, needed to 
validate the thermodynamic model are shown in Section 
4. The bubble growth continues even when the probe is 
moving outwards from the Sun (> 2.8 AU, grey area). Ob¬ 
viously, at larger distances the bubble growth slows down, 
since the probe is being bombarded by the smaller proton 
fluxes, see Eq. 2. 


For simplicity it is assumed that a single bubble is a 
half of a sphere with a radius of rp The gas within a 
bubble behaves to a good approximation like an ideal gas: 

N 

p i ^ = £jV Ha ,ijfcBT, (5) 

j 

where p; is the pressure of the gas, kp is the Boltzmann 
constant, T denotes the temperature of the sample, N is 
the number of time steps up to a given state of bubble 
growth, hence the irradiation time of the sample after N 
time steps is N x Afj. 

The number of recombined H atoms, subtracted by 
those which diffuse from the sample out (iVdiffj) is: 

Ahj = / E AijA(l — BS) + iVdiffj, (6) 

V diff ,j = -D H (T)-^f-(A-A h ) Afj. 

The term Dn(T) determines how much Hydro¬ 

gen diffuses from the sample out (per unit area and time). 
Therefore, by dimension it is a flux of outflowing Hydro¬ 
gen from the specimen. Its constant value results from the 
following fact. The model assumes that the sample is ex¬ 
posed to the protons having kinetic energies from eVs to 
8 keV. According to the data (see, Fig. 1), the magnitude 
of the proton flux remains almost constant for the con¬ 
sidered energy range. Therefore, the sample is uniformly 
populated by the protons (H atoms after recombination) 
to a depth of dpR. The dpR is calculated for the fastest 
protons. Therefore, there is no differentiation of Hydrogen 
concentration with respect to the depth. 

In Eq. 6 A is the area of the sample irradiated by 
the protons, A b is the area of the sample covered by the 
bubbles. BS is the factor of backscattered ions. If BS is 1 
then all of the incident ions are backscattered. If BS is 0 
then all of the incident ions penetrate the target. Dp(T) 
is the diffusion coefficient for H atoms in a given material, 
(qpj is the number density of H atoms which may diffuse 
through the lattice in the j th period of time, dpp(E) is the 
so-called projected range. It is defined as an average value 
of the depth to which a charged particle will penetrate in 
the course of slowing down to rest. This depth is measured 
along the initial direction of the particle, and it depends 
on the kinetic energy of the particle (Berger et ah, 2005). 

The number of Hydrogen molecules added in the j th 
period of time to the i th bubble is constant and 

given by: 


3.2. Thermodynamic approach to blistering process 

In the following, a thermodynamic model of bubble 
growth is proposed. The model is based on the assump¬ 
tion that the growth proceeds quasi-static i.e. during a 
jth p er i oc i 0 f time Afj a small portion of H 2 -molecules, 
7V H2 .i,j, is added to the i th bubble and a thermodynamic 
equilibrium is rapidly re-established. 


iVka.iJ = 0.5Gi-/VHj f)max(s) (7) 

53 G i = !, dimG = Ng. 

i 

Here 0.5 denotes that a single H 2 molecule consists of two 
H atoms, G is a matrix, its role is to redistribute cer¬ 
tain number of H 2 molecules into the bubbles. The input 
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pattern of H 2 molecules into the bubbles follows experi¬ 
mental findings of Kamada et al. (Kamada et al., 1986). 
IVg is the total number of bubbles on the irradiated sam¬ 
ple. While 100% of protons recombine into H atoms in 
the metal lattice, only a part of them recombine to H 2 
molecules (Canham et al., 1989). Hence the ? 7 max (s) co¬ 
efficient is the ratio between the number of H 2 -molecules 
and the H-atoms in the lattice. The H 2 molecule is formed 
when electrons of two H atoms have anti-parallel spin s, 
otherwise the molecule cannot be created. Therefore, at 
most half of the H atoms can form H 2 molecules, hence 
i)max(s) = 0.5. Not all of the H 2 -molecules will merge into 
H 2 -clusters and finally form H 2 -bubbles. Thus, the coef¬ 
ficient £ denotes the ratio of the number of H 2 -molecules 
inside and outside the bubbles. 

The first step to estimate the radius of the i th bub¬ 
ble is to calculate the Helmholtz free energy of the whole 
configuration, F con fi g . Since the free energy is an additive 
quantity, the total free energy of bubble formation is a 
sum of following quantities: free energy of H 2 gas inside 
the i th bubble (F gaS) i), of the metal surface deformation 
(F m d.i) caused by the bubble growth itself, of the surface 
free energy (F sur f^) of the bubble cap, of the free energy 
of H 2 -molecules (Fh 2 ) and of H-atoms (Fh) placed outside 
the bubbles but within the metal lattice. The Helmholtz 
free energy of the whole configuration described above is 
then: 


N b 

Fconfig = E (^& as l + ^md,i + ^surf.i) + F Ha + F H . (8) 
i 

The next step is to estimate the free energy of the i th 
bubble. It consists of the free energy of the gas filled in 
the bubble, the free energy of metal deformation, and of 
the bubble cap surface free energy. 

Using the thermodynamic relation between gas pres¬ 
sure and its Helmholtz free energy p = (jp/) T together 
with the equation of state Eq. 5, the free energy of a gas 
within the i th bubble is: 

F gas ,i = - £ Nh 2 ..ifeTln (, (9) 

\ Mnin J 
J 

where U maXj i is the maximum volume of a given bubble. 
The model assumes that two H 2 molecules form the small¬ 
est (’’initial”) possible bubble, its volume is denoted by 
U m in- The radius of such a bubble is approximately 1.45 
Bohr radii (Ree & Bender, 1979). Every bubble will crack 
if the pressure of the gas inside is higher than the pressure 
exerted by the metal deformation of the cap. The relation 
between the pressure of the gas, surface tension a, and the 
bubble radius corresponding to Vmax.i is (Lautrup, 2011): 

2(7 

P gas, insite bubble Poutside bubble — — • V-LUJ 

^max,i 


Since the sample is placed in vacuum, the pressure outside 
the bubble is set to zero. 

The free energy of metal deformation F m( j,i caused by 
the gas pressure inside the bubble with radius r\ can be 
found in (Landau & Lifshitz, 2009), and is given by: 


F m d,i 


47r rf(l + 7) 2 

3 E y Pi 


( 11 ) 


Here 7 is the Poisson coefficient, i.e. ratio of transverse 
to axial strain of a sample material, Ey is the Young’s 
module. 

The free energy of a surface of a cap of the i th bubble 
is given by (Martynenko, 1979): 


F surf , i = 47rrf <7(T). (12) 

The Helmholtz free energy of the H 2 -molecules located 
at certain positions in the metal lattice but outside the 
bubbles is calculated in the form F = F; nt — TS. Where 
Fi nt is the internal energy of molecules/atoms located at 
certain positions in the metal lattice. Applying the statis¬ 
tical definition of the entropy S, this free energy is: 


/ N \ 

Fh- 2 = ( (I 3 ) 

N 0 ) J ’ 

where N^ 2 is the total number of H 2 molecules inside the 
sample, eu 2 is the binding energy of H 2 molecule to a va¬ 
cancy. A detailed derivation of the Eq. 13 is presented in 
Appendix A.l. Nq is the number of lattice sites, which 
can be expressed by: 

No = -/Va^pr——, (14) 

1 Vl u 

where Na is the Avogadro’s number. M u is the molar 
mass of the sample’s material. 

The Helmholtz free energy of H atoms located at cer¬ 
tain positions within the metal lattice is: 


£h 2 + k B T I 11 


F h = (N* - 2 Nl) 


£h + fcsFln 


N£ - 2 Nl 

No 


_ (15) 

where en is the migration energy of the H atom in the 
metal lattice, and N ^ is the total number of H atoms in 
the sample. A detailed derivation of the Eq. 15 can be 
found in Appendix A.2. 

Since now each term of Eq. 8 is determined, the next 
step is to estimate the radius r, of the i th bubble at given 
time t. This will be achieved by assuming that the process 
of bubble growth proceeds quasi-static in thermodynamic 
equilibrium: 
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dF c 


config,i 




= 0. 


(16) 


This condition leads to the following fifth order equation 
for Ti: 


87rS; 


:i a(T)rf - H ir f + ^ ( E f JV Ha , U ) 
X [2N n - 33;,j Ef ^H 2 ,i,j] = 0, 


k%T 2 


(17) 


H;j is defined below in Eq. 19, Hi denotes the abbrevia¬ 
tion: 


d-Fgas,i _ dFii _ dFg 2 

dNK 2 ,i,j dNft 2 ,i,j dN H2<i j 


(18) 


The derivatives of the Helmholtz free energy, of the gas 
inside the i th bubble, of metal deformation caused by the 
bubble, of the i th bubble cap surface, of H 2 molecules, and 
of H atoms with respect to the number of H 2 molecules 
added at the j th time step to the i th bubble are presented 
in Appendix A.3. 

A realistic model of bubble radius growth, can 

be estimated by following Gedankenexperiment. Obvi¬ 
ously at the beginning of the bubble growth process, the 
differential increase of the bubble radius is higher than at 
its end. It is implied, that the number of H 2 molecules 
in the system is conserved and at each time step one of 
them merge into a bubble. After At the bubble consists of 
2 H 2 molecules, hence the number of molecules increases 
by 50%. At the time 2Af the bubble consists of 3 H 2 
molecules, hence the number increase is now 33.3%, and 
so on. Therefore S is: 


the foil has not been exposed to the electromagnetic ra¬ 
diation and/or charged particles, and is considered to be 
a perfect mirror with the reflectivity of R = 1. It means 
that all of the incident light rays are reflected perfectly, no 
light ray is absorbed or diffusively reflected by the target. 
Later, when the foil has been irradiated by a flux of pro¬ 
tons and molecular Hydrogen bubbles have been formed 
on its surface, the reflectivity of the degraded foil will be 
reduced. This deterioration is calculated in the following 
way: the foil is covered by a grid with a fixed single cell 
size of e ce ii x e ce ii, see Fig. 4. The reflectivity of a single cell 
is by definition , where A q is momentum transfer of 

a photon to the i th cell of the degraded foil, while Ag maXj ; 
is the momentum transfer of a photon to the i th cell of a 
perfect mirror. 



Figure 4: A fraction of the foil with one spherical bubble is shown. 
The size of a single cell of the grid is e ce n x e ce ii. 


A n 


A Ni 




= 3 n.o, 



(19) 


The exponent a is a model parameter of the bubble growth. 
The value ^ corresponds to the Gedankenexperiment pre¬ 
sented above. However, the true value of the a parameter 
differs from that. In the process of bubble growth, par¬ 
ticles (the Hydrogen) are added to the system i.e. the 
probe is permanently irradiated by the protons, they pen¬ 
etrate the target and recombine to the Hydrogen. On the 
other hand, both, due to the diffusion process and bubble 
cracking, some Hydrogen atoms leave the system. There¬ 
fore, the number of Hydrogen atoms in the system is not 
conserved. Hence, a series of experiments have been per¬ 
formed to estimate a realistic a parameter, results are pre¬ 
sented in Section 4. 


3.3. Reflectivity of a metallic foil covered with bubbles 
The momentum transfer of a photon to an ideal reflect¬ 
ing surface is given by A q = 2qcos9, where the factor 2 
is just in accordance with specular reflectivity. Certainly, 
the surface quality will suffer during the irradiation with 
protons from progressing bubble formation. At time t = 0 


Therefore, taking into account all cells, one has: 


A R = 


Ef ^ 11 Ag; 

Ef- 1 Ag ma: 


( 20 ) 


Here jV ce n is the number of cells. The path of photons is 
directed parallel to the foil surface normal. Therefore, at 
time t = 0 the foil was a perfect mirror without surface 
imperfections and 9{ = 0. Later, when the surface is popu¬ 
lated with bubbles, 9\ will vary between 0° and 90°. Thus, 
Eq. 20 reduces to: 


A R = 


Er- 2? cos 9{ 
iVcell X 2 q 


Ef- cos 9j 
Acell 


( 21 ) 


4. Results 

A following set of experiments were performed. Three 
probes (Al, A2, and A3) were exposed to a flux of 2.5 keV 
protons, each one with longer irradiation time, see Table 
2, where ig is a number of days in space until a probe will 
collect a given dose of protons. Results are shown in Fig. 5. 
From top to bottom, the pictures correspond to the probes 
Al, A2, and A3, respectively. Average sizes of bubbles 
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have been estimated to 0.18 ± 0.05 /im, 0.19 ± 0.05 /mi, 
and 0.2 ±0.05 /im for probe Al, A2, and A3, respectively. 
There is a strict correlation between a dose of protons and 
the average bubble size for a given population. The higher 
the proton dose, the larger the bubble sizes. Examining 
the electron microscope pictures, the surface density of 
bubbles has been estimated to ~ 10 8 cm -2 . 

We have confirmed formation of molecular Hydrogen 
bubbles by a number of experiments. We have varied pro¬ 
tons flux, energy, and specimens’ temperature. Bubbles 
were observed for specific parameters range. Fig. 3 repre¬ 
sents our findings. 

The here used three experimental findings are the mini¬ 
mum number of specimens necessary to validate the model. 


Table 2: Test parameters 


Probe symbol 

T[K] 

E [keV] 

D [p + cm ' z ] 

is [days] 

Al 

323.0 

2.5 

7.8 x 10 17 

4.8 

A2 

323.0 

2.5 

8.2 x 10 17 

5.0 

A3 

323.0 

2.5 

1.3 x 10 18 

7.9 


4-1- Bubble growth law 

Bubble growth in blistering process is dominated by 
diffusion of hydrogen atoms into the bubbles. A simi¬ 
lar problem of the bubble growth (Ostwald ripening) was 
widely discussed in literature (e.g. Chaikin & Lubensky, 
2000; Jones, 2002; Lifshitz & Slyozov, 1959; Lifshitz, 1962; 
Rutenberg & Bray, 1995; Wagner, 1961). It has been found 
the power law f(t) ~ ts for the average bubble radius 
f in the systems with the total mass being a conserved 
quantity. This type of the systems exhibits diffusion lim¬ 
ited growth and is well described in terms of the Lifshitz- 
Slyozov-Wagner (LSW) theory (Lifshitz & Slyozov, 1959; 
Wagner, 1961). The Hydrogen bubbles follow the diffusion 
limited growth and the question can be raised whether the 
LSW theory applies to them. In this case the growth rate 
of the bubble average radius f(t) could be represented by 
the following phenomenological equation: 


dr(t) = £b_ L, , _ ac 
dt f(t) [ [> f(i) 


( 22 ) 


similarly as in (Lifshitz & Slyozov, 1959), where D B is the 
diffusion coefficient of Hydrogen atoms, 


ac 


2<jV ia C 00 

k B T 


(23) 


is a parameter which depends on the surface tension a of 
a bubble, the Hydrogen atom volume V mi the concentra¬ 
tion Coo of Hydrogen atoms far beyond the bubble and 
temperature T, A (t) is the difference of concentration of 
Hydrogen atoms at the bubble boundary and Coo- In the 
LSW approach, the total mass of a system is a conserved 
quantity, in consequence ' — const, and the solution 

of Eq. 22 takes the form of f(t) ~ t 1 / 3 . 



S pm 




Figure 5: Electron microscope pictures of probes Al (top), A2 (mid¬ 
dle), and A3 (bottom). 















The three probes Hi, A 2 , A 3 from the pictures in Fig. 
5 were exposed to a constant proton flux. Taking into ac¬ 
count that the amount of the Hydrogen atoms produced 
by the proton flux is related directly to the existing elec¬ 
tronic structure of the irradiated probes as well as that 
the thickness of the superficial layer including the molecu¬ 
lar Hydrogen bubbles is relatively constant in size (it only 
shifts deeper into the sample after the superficial bubbles 
break) then the system of the growing bubbles resembles 
the total mass conservation system. Hence, it could be ex¬ 
pected that at least the largest bubbles which appear on 
the irradiated probes surface follow the growth law t 1 / 3 as 
well as they undergo coalescence phenomena. The latter 
property is evident, e.g. for probe A 3 in Fig. 6. Concern¬ 
ing the question of the possible existence of the t 1 / 3 scaling, 
the average bubble radius, respectively for the probe Hi, 
A 2 , A 3 (Table 2) has been divided by the duration of the 
proton flux exposure in the power of 1/3, i.e., 4.8 1 / 3 , 5 1 / 3 
and 7.9 1 / 3 . All scaled radii take the same value close to 0.1. 
In addition, in Fig. 6, the distribution of the bubble radii, 
r/t 1//3 of the growing bubbles has been plotted. It can be 
observed that the largest bubbles seem to follow the f 1 / 3 
power growth law because their scaled distributions coin¬ 
cide. The left hand side of the bubble radius distribution 
is influenced by the bubbles which were formed at later 
time moments and it is the reason for the strong deviation 
of their size from the t 1 / 3 power law. It can be concluded 
from the surface analysis of the probes Hi, A 2 and A 3 that 
their degradation after they were exposed to proton flux 
is not faster than it is predicted by LSW theory. However, 
it is necessary further investigation on the problem both 
experimentally and theoretically. 

4-2. Validation of the thermodynamic model 

For numerical simulation a 10 fxm x 10/im foil was spec¬ 
ified. That choice allows to simulate a smaller number 
of bubbles, i.e. it decreases the computation time of the 
simulation. It implies also an important assumption that 
surface arrangement of the bubbles is isotropic i.e. any 
10 /im x 10/iTO area of the irradiated sample is indistin¬ 
guishable. Table 3 collects all of the model parameters 
used in the simulation. The first set of parameters charac¬ 
terize mechanical and thermo-optical properties of vacuum 
deposited Aluminum on UBS’s Upilex — S® foil. Second 
set specifies values of the parameters which have been used 
to fit the model to the experimental data presented here. 

To fit the proper gradient of bubble growth, the a pa¬ 
rameter was set to 0.6, see Eq. 19. Comparison of the aver¬ 
age bubble size of the experimental and numerical findings 
are drawn in the top plot of the Fig. 7. The £ parameter 
was set to 0.98. It determines the height of the curve. A 
decrease of the specular reflectivity of the foil as a function 
of time is shown in the bottom plot of the Fig. 7. The 
decrease of the reflectivity is 3.0, 3.2, and 4.6 % in compar¬ 
ison to the non-irradiated foil for t s = 4.75, 5.0, and 7.9 
days, respectively. Clearly, the larger the bubble sizes, the 
larger the specular reflectivity decrease A R in comparison 



bubble radius, r [ji m] 



1/3 

scaled bubble radius, r/t [\i m] 


Figure 6: Top plot, distribution of bubble radius r for bubbles in Fig. 
5, respectively for probe Al, A2 and A3. Bottom plot, the same as 
top plot but for the scaled radius, r/t 1 / 3 . 
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to the non-irradiated foil. A distribution of the bubbles at 
three different time steps: is = 4.8, 5.0, and 7.9 days is 
shown in Fig. 8. During time, the probe collects higher 
dose of protons, therefore, the distribution drifts i.e. size 
of bubbles increases. At time is = 7.9 days most of the 
bubbles have sizes in range 0.19 to 0.24 pm, there are only 
a few which have sizes larger than 0.25 pm. 


Table 3: Model parameters 


Symbol 

Value 

Q 

2.7 [g cm -3 ] 

M u 

26.98 [g moH 1 ] 

Ey 

69 x 10 10 [dyn cm -2 ] 

7 

0.33 

£H 

0.52 [eV] (Linderoth, 1988) 

£ h 2 

0.06 [eV] (Lu & ICaxiras, 2005) 

as 

0.093 

£t 

0.017 

BS 

0.02 (Ziegler, 2013) 

A 

100 [/cm 2 ] 

T 

323 [K] 

'Hma x(^) 

0.5 


0.98 

a 

0.6 

n b 

10 8 [cm -2 ] 




Figure 7: Time evolution of an average bubble radius from the popu¬ 
lation (top plot), specular reflectivity decrease due to bubble growth 
(bottom plot). 


5. Conclusions 

It has been proven that thermodynamic model is a flex¬ 
ible tool to simulate and to reproduce the real growth of 
the molecular Hydrogen bubbles. However, the estimated 
a and £ parameters are applicable only for the here pre¬ 
sented experimental findings. These parameters depend on 
type and temperature of the irradiated material. There¬ 
fore, change of the material type and the experimental 
conditions requires an update of the model parameters. 

The time evolution of decrease of the specular reflec¬ 
tivity A R is a model prediction. It is highly required to 
perform experimental confirmation of that findings, since 
the real reflectivity decrease can differ from that estimated 
here. 

For sure, the here presented experimental findings show 
that physical parameters of the metal surface, e.g. rough¬ 
ness, reflectivity, are significantly changed under condi¬ 
tions as prevalent in the interplanetary space. There¬ 
fore, statement present in the ECSS standard (ECSS-Q- 
70-71A, 2004, p. 51) that metals do not suffer from space- 
environment is not longer valid. 

The thermodynamic model requires further improve¬ 
ments. The considered aging factor, the solar protons, are 
not the only one which can influence the bubble growth 
process. The solar wind is also essentially made up of elec¬ 
trons and small proportion of heavier ions (Meyer, 2007). 
Additionally, electromagnetic radiation in the lower wave¬ 
length range has also to be taken into consideration. These 


£20 


Al, t=4.8 days 
A2, t=5.0 days 
A3, t=7.9 days 


0.20 0.22 0.24 

bubble radius r [/im] 


50 

JS 

S 40 

D 

-Q 

o30 

ai 

-Q 

E20 

C 

10 


Al, t=4.8 days x 
A2, t=5.0 days • 
A3, t=7.9days • 



0.09 0.10 0.11 0.12 0.13 0.14 0.15 

scaled bubble radius r/t 1/3 [^m] 


Figure 8: Bubble size distribution at a 100 /Ltm 2 sample at three 
different time steps: 4.8, 5.0, and 7.9 £g. 
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degradation factors can slow down the bubble growth. The 
growth deceleration can be explained as follows. Hydrogen 
molecule dissociate at the energy of 4.5 eV (Balakrishnan 
et ah, 1992). The dissociation may be caused by the UV- 
light with the wavelengths < 274 nrn. The H 2 gas within 
the bubbles can then be partially dissociated, and H atoms 
can diffuse easily through the bubble caps. As a result the 
bubble growth process may slow down. The deceleration 
can be strengthen by heavier ions generated by the Sun 
e.g. a-particles. Their diameter is much larger than that 
of protons or electrons, hence, collisions between the H 2 
molecules and the a-particles within the bubbles can ad¬ 
ditionally increase the dissociation efficiency. 

The present condition for the bubble crack mechanism, 
Eq. 10, assumes that the pressure outside the bubbles 
is negligible small. Under the real space conditions the 
electromagnetic radiation will exert a pressure on the caps, 
hence, their sizes may be smaller. On the other hand, 
bubble caps loose thermal contact with the base material 
and they become overheated (Astrelin et al., 2010). As 
a consequence the caps can brake and the H 2 gas can be 
released. That aspect of the blistering process needs to be 
examined. 

By these reasons further experimental studies are planned, 
e.g. the mentioned influence of the UV-light on the bubble 
growth dynamics. 

6. Acknowledgements 


of H 2 molecules in the i th bubble. The second summa¬ 
tion over number of bubble i th counts the number of H 2 
molecules within all bubbles. Therefore, performing the 
subtraction one gets the total number of H 2 molecules out¬ 
side all bubbles, located at certain positions in the metal 
lattice. The entropy k B In f 1 of the collection of the H 2 
molecules placed on a lattice sites is (Seitz, 1940): 


(*S a - E?* Ej N *h 2 ..j) ! [no - - E^ Ef JVHa.u)] ! ’ 

^h 2 - Ef® Ef 

N 0 

(A.2) 

The internal energy E- m t is given by the following relation: 

A3 N \ 

> ( A -3) 

where the en 2 is the binding energy of the H 2 molecule to a 
vacancy (Lu & Kaxiras, 2005). The internal energy of H 2 
molecules located in the metal lattice sites is a product of 
the binding energy of a single H 2 molecule and the number 
of molecules. 

The Helmholtz free energy of the H 2 molecules located 
outside the bubbles at certain positions in the metal lattice 
is then: 


< 2 -EE%.i- j 


In 
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Appendix A. Bubble formation 


N 

F H 2 = I N H 2 ~ 

i j 


(A.4) 


£h 2 + /cbUIii 


An 


Appendix A.l. Helmholtz free energy of H 2 molecules placed 
in certain positions in the metal lattice 
To get the relation of the Helmholtz free energy of the 
H 2 molecules located at certain positions in the sample 
but outside the bubbles, the definition of the free energy 
is used: 


Appendix A.2. Helmholtz free energy of H atoms in the 
sample 

The number of H atoms in the sample and outside the 
bubbles is: 


F = E int - TS , S = k B Inn, (A.l) 

here E- lnt is the internal energy of the H 2 molecules located 
at certain positions in the metal lattice, S denotes the 
entropy, while Q represents the number of ways in which 
the H 2 molecules can be arranged on the Nq lattice sites 
(Eq. 14). 

The number of H 2 molecules located outside the bub¬ 
bles at certain positions in the metal lattice is A ^ 2 — 

Ef Where A ^ 2 is the total number of H 2 

molecules in the sample, while ^E® Ej V AH 2 ,i,j is the total 
number of H 2 molecules within all bubbles. The summa¬ 
tion over the number of time steps j counts the number 


n b n \ 

EE^ij + Agf bubb H, (a. 5 ) 

where Ag is the total number of H atoms in the sample, 
so the number counts all of the incident Hydrogen ions 
which have recombined into Hydrogen atoms. Some of 
the Hydrogen atoms have recombined to H 2 molecules and 
some of the molecules are forming the bubbles. Hence to 
get the number of H atoms located on the lattice sites one 
has to subtract the total number of Hydrogen atoms 
and those Hydrogen atoms which build H 2 clusters and 
H 2 bubbles. The reason of the factor 2 is that a single H 2 
molecule consists of two H atoms. 
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The procedure to estimate the Helmholtz free energy 
of H atoms in the sample is the same as in Eq. A.4. Hence 
the term is: 


^-^surf ,i 


cW H2 


= 8 lrn f)N 1 a T)- 


(A.ll) 


E h = N B + k B T In > (A.6) 

where e B is the migration energy of the H atom in the 
metal lattice. The migration energy is defined as the min¬ 
imum energy which has to be added to the H atom in order 
to remove it from the lattice site. 

Appendix A.3. The derivatives of Helmholtz free energy 
of: gas of the i th bubble, metal deforma¬ 
tion caused by the bubble, surface of the 
i th bubble cap, H 2 molecules, and H atoms 
located outside the bubbles. 

The equlibrium condition of the process of the i th bub¬ 
ble growth is: 


dF c 


config,i 


5ATH 2 ,i,j 


= 0 , 


(A.7) 


The derivative of the Helmholtz free energy of the H 2 
molecules (Eq. A.4) located outside the bubbles at cer¬ 
tain positions in the metal lattice is: 


OF h „ 


OATh, 


= -iViVg . 

£h 2 + ksT 

1 +ln 

N 

{ 





- Ei B £f Nh 2 


N 0 


(A. 12) 

The derivative of the Helmholtz free energy of the H atoms 
(Eqs. A.6) located at certain positions in the metal lattice 
is: 


dF H 

9An 2 ,i,j 


= —2 M b t 

N 

/ NS - 2NZ Y 


(A.13) 
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The assumption is fulfilled when the time scale of the 
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The free energy of a metal deformation caused by expand¬ 
ing i th bubble is given by the Eq. 11, hence its derivative 
is: 
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